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Abstract 

We study the NGS (Non-dissipative Gravitational Singularity) model, which 
successfully describes the non-linear stage of evolution of perturbations (see 
[1], [2] and references therein). This model predicts DM density distribution 
p[r) ~ r~ a with a ~ 1.8 which holds from very small distances r m i n ~ 0.01 pc 
up to very large distances r max ~ 5 Mpc. Assuming the neutralino to be a 
CDM particle, we calculate the annihilation of neutralinos in the vicinity of 
the singularity (Galactic Center). If neutralinos are the dominant component 
of DM in our Galaxy, the produced energy is enough to provide the whole 
observed activity of the GC. Neutralinos of the most general composition 
and of mass in the range 20 GeV < m x < 1 TeV are considered. We 
find the neutralino compositions which give the relic density needed for the 
Mixed Dark Matter (MDM) model and we evaluate for these compositions 
the high-energy (i? 7 > 100 MeV) gamma-ray flux under the constraint that 
the radio flux is lower than the observational limit. The compositions with 
the detectable gamma-ray flux which we found are provided by a set of 
almost pure gaugino states with the neutralino mass between 100 and 500 
GeV. 

We demonstrate that a detectable high-energy gamma-ray flux is pro- 
duced by the neutralino annihilation also in the case when neutralinos provide 
a small fraction (down to 0.1%) of the DM in our Galaxy. The predicted flux 
is ~ 10" 7 - 10" 8 cm" 2 s" 1 for £L > 300 MeV 



- 1 - 



1. Introduction. 

The spectrum of perturbations as observed by COBE, IRAS and CfA is 
most naturally explained by the Mixed Dark Matter (MDM) model [3] with 
~ 60 - 70% of Cold Dark Matter (CDM), ~ 20 - 30% of Hot Dark Matter 
(HDM) and ~ 10% of baryonic matter . The most plausible candidates for 
the MDM particles are the r-neutrino as a HDM particle and the neutralino 
X as a CDM particle. According to this model 5p/p on the galactic scale 
A ~ 1 Mpc is dominated by the CDM particles, while on the scales A ~ 
10 — 100 Mpc 5p/p is dominated by the HDM particles. The baryons, which 
decoupled from the radiation much later than CDM particles, streamed into 
the potential wells built by the space distribution of CDM particles. As 
a rough estimate of the average fraction of the baryon density in the halo 
one can take the relative density of baryons in the whole Universe (~ 10%) 
and this assumption does not contradict the recent observations of MACHO 
candidates [4]. 

When Sp/p on the galactic scale approaches one, the evolution of fluc- 
tuations enters the non-linear stage. There are two solutions for this stage, 
determined by the space symmetry of the initial fluctuations, i.e. of those 
with 5p/p ~ 1. The first solution was found by Zeldovich [5] for a case of 
asymmetric initial distribution of matter relative to the local maximum of 
density. In this case the evolution terminates with the production of the 
flat formations - the Zeldovich pancakes. For this solution the potential well 
within a galaxy is not formed and the CDM gas is not self-trapped. 

The other solution, found first in ref. [1], exists if the initial distribu- 
tion is more symmetric. In this strong singularity is formed and the 
CDM gas is self-trapped. Each singularity, which is called Non-dissipative 
Gravitational Singularity (NGS), is associated with a galaxy. The main pre- 
dictions of this model according to ref. [2] , can be summarized as follows: 

i) The density distribution of CDM particles centered at a NGS is given 

by 

p(r) oc r~ a (1) 

with a ~ 1.8. The distribution (1) is valid down to very small distances 
r min where a cut off can be caused by several phenomena[2]. The largest 
value of r min is related to the case of a black hole in the NGS. For a black 
hole with a mass M ~ 10 6 M this value is r min ~ 0.01 pc. The maximum 
value up to which the distribution (1) holds is not predicted by the NGS 
model, because this is the initial condition for the non-linear problem. We 
suggest that the value of r max ~ 5 Mpc, as it is given by the two-point 
correlation function, is determined by the scale where HDM dominates in 
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the perturbation amplitudes. The distribution (1) is in excellent agreement 
with the observations. 

ii) The predicted density of CDM at the distance Tq = 8.5 kpc from the 
Galactic Center is equal to Pcdm — 0-2 GeV, if we take < pcdm >~ 0.7 p c 
(p c = 1.88-10~ 29 /i 2 g cm -3 ). This result is in a good agreement with the DM 
density obtained from the rotational curves in our Galaxy, 0.3 GeV cm -3 < 
Pdm < °- 43 GeV cm~ 3 : for the total DM density. 

iii) Annihilation of CDM particles in the vicinity of the singularity results 
in a pointlike radiation source in the Galactic Center. 

In ref.[2] the luminosity due to the neutralino annihilation in the NGS 
model has been evaluated, together with the ensuing gamma-ray and radio 
fluxes. Comparisons of the calculated fluxes with the available experimental 
upper limits were then used in ref.[2] to obtain lower limits for neutralino 
masses in the case of almost pure neutralino configurations (either gauginos 
or higgsinos) and under the hypothesis that neutralinos provide a fraction 
of about 70% of DM density in the Galactic halo. In ref.[2] the theoretical 
evaluations referred mainly to not too heavy neutralinos, i.e. m x < 0(mw)- 

In this paper we present a more general analysis valid for a wide range of 
neutralino masses ( 20 GeV < m x < 1 TeV) and of neutralino compositions 
and we mainly address the question: can the annihilation of neutralinos 
around the singularity result in detectable gamma-ray flux at energy E 1 > 
100 MeV? Apart from the favourite case when the neutralino is the only 
CDM particle, we shall study the case of arbitrary Q x = p x j p c , which in the 
framework of the MDM model implies the presence of other CDM particles 
X with O x > O x 

2. Neutralino 

We shall consider the neutralino in the context of the Minimal Super- 
symmetric Standard Model (MSSM) . The neutralino is defined as the lowest- 
mass linear combination 

X = ZuW 3 + Z 12 B + Z 13 H U + Z 1A H d (2) 

where W% and B are the SU(2) and U(l) neutral gauginos and H U} d are the 
higgsinos. The space of neutralino states is determined by three independent 
parameters: the masses of B and W 3 , Mi and M 2 , respectively, with the 
GUT relation between them M\ ~ 0.5 M 2 , the Higgs mixing parameter p,, 
and the ratio of the two v.e.v.'s which give masses to up-type and down-type 
quarks, tan/3 = v u /vd- The natural range for tan/3 is: 1 < tan/3 < m t /mb- 
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In the calculations discussed below the parameters M 2 and /i are varied in 
the wide ranges: < M 2 < 6 TeV and -3 TeV < /x < 3 TeV. 

The neutralino is assumed to be the lightest SuSy particle and to be 
stable because of R-parity conservation. The decoupling of neutralinos in 
the early Universe with the subsequent annihilation results in the neutralino 
relic abundance 

^ = " !1 - W - 11 ®)'(2^)'<(S^)- < 3 » 

where h is the dimensionless Hubble parameter, Xf = Tf/m x , Tf is the neu- 
tralino freeze-out temperature, T 1 and T x are the present temperatures for 
relic photons and neutralinos, Nf denotes the effective number of degrees of 
freedom at the freeze-out temperature. Nf ~ 100 for m x > m\y- The ther- 
mal average of the xx-annihilation cross-section < av > at the temperature 
T = xm x is parameterized in (3) as 

< av >= a + bx . (4) 

The expression ax f + 1/2 bxj- =< av >i nt in eq.(3) represents the integrated 
cross section from Tf to the present temperature. 

Since we wish to discuss neutralinos with masses up to 1 TeV, many 
channels have to be considered in the neutralino annihilation processes. In 
the present paper, for the values of < av > and of the neutralino relic abun- 
dance, we employ the results of the thorough analysis of ref.[6], where the full 
set of possible final states have been considered: fermion-antifermion pair, 
pair of two (neutral and charged) Higgses, one Higgs boson-one gauge boson 
states and pair of two gauge bosons. In the evaluation of a, the exchanges 
of the following particles : Z boson, Higgses, sfermions /, neutralinos and 
charginos, have been included depending on the specific final state. As for the 
values of the masses of the neutral Higgs bosons, the two CP-even: h, H (of 
masses rrih and ran-, with rrih < mg) and the CP-odd: A (of mass ttia), we 
employed the usual relationships which include radiative corrections. These 
imply that if, for instance, rrih is used as a free parameter, then rriA and 
run turn out to be dependent on rrih, tan/? and also (through the radiative 
terms) on rrit (top-quark mass) and rh (mass of the top scalar partners, con- 
sidered here as degenerate). We remind that the LEP lower limit for rrih is 
m h > 50 GeV [7]. 
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In our further calculations we give most attention to those compositions 
of x which provide fi x ~ 0.6 — 0.7. However, we shall also consider of 
smaller fi x , which implies that neutralino compose only a fraction of CDM 
in the Universe. 

3. The pointlike source in the Galactic Center due to xx-annihilation 

For the distribution (1) annihilation between neutralino pairs would be 
very efficient in a region of linear size ~ r m i n ~ 0.01 pc which implies for 
the source the angular size ~ 10~ 6 — 10~ 7 rad. The total luminosity of the 
source is 

L = 2m x < av > P (5) 

where P is 



^4 ( * c)2r| (^r 6 <6> 

and < av >o denotes the value of < av > at the present temperature. 

Here we have employed for the neutralino density profile in our Galaxy 
the NGS distribution of Eq.(l) and have denoted by p l ° c the local (solar 
neighborhood) neutralino density. Due to different space distribution of 
CDM and HDM in the Universe, we expect p^/Pdm ^x/^- However, 
for simplicity we used in our calculations P x c /Pdm = ^x/^- For r min we 
shall use the capture radius of the black hole with mass M ~ 10 6 M© [2]: 
r min ~ 0.01 pc. 

In the evaluation of cross section < av >o in Eq.(5) we use the ap- 
proximation < av >o— a, since at present, inside the Galaxy, neutralinos 
most probably have a gaussian velocity distribution with an average value 
v ~ 300 km s _1 and consequently x ~ v 2 /6 ~ 10 -7 . Annihilation of neutrali- 
nos into quarks and leptons thorough the main channels 66, tt, W + W~, and 
ZZ has been taken into account. The part of luminosity due to fermion pro- 
duction is dominated by heavy (b and t) quarks. The 6(6) quarks, either pro- 
duced directly or by t(t) decays, do not decay within the confinement region, 
but produce copiously lighter quarks which then turn into pions and kaons. 
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For the calculation of the hadronization process we used the Lund Jetset 
7.2 Monte Carlo simulation program [8]. This program, for any given initial 
configuration (i.e., the x~ X annihilation final state F), allows one to follow 
the hadronization and to infer the electron and 7 final distributions. We have 
performed this analysis for each final state mentioned above. The statistical 
sample was taken to be 2 • 10 4 events and the relevant physical parameters for 
each process were taken from the PDG [9]. It is worth mentioning that the 
photon and electron spectra computed at different center-of-mass energies 
E crn ~ 2m x , when plotted in terms of the variable x = E /(E cul /2), exhibit 
remarkable scaling properties[10]. Thus xx-annihilation in the GC results 
in high energy (E 1 > 100 MeV) gamma-radiation, X-rays (mainly from the 
Compton scattering of the electrons), in gamma-ray lines (from XX ~ > 77 
and e + e~ — > 77 annihilation) and in radio-emission (from synchrotron ra- 
diation of electrons in magnetic fields). The estimates in ref.[2] show that 
expected high energy (E 7 > 100 MeV) gamma-radiation and radio-emission 
can be at the level of the observable fluxes or may exceed them. 

The quantity of interest for observations is the integral gamma-ray flux 



where P is given by (6), x = E^/m x , < av >f is the XX annihilation cross 
section in the final state F and dN 1 {x)/dx is the photon energy spectrum 
mainly due to neutral pions decay. The existing limit [11] on gamma-ray 
flux from the GC is related to E y ~ 300 MeV. More generally, for future 
experiments E 1 of interest is ~ 70 — 100 MeV. 

The energy flux for radio-emission at frequency v at Earth can be cal- 
culated as 



where E e is the electron energy at generation, E is the electron energy at 
the moment of synchrotron radiation, b(E) = (dE/dt) s is the synchrotron 
energy-loss of an electron and P{E, v) is the energy-loss of electron with 
energy E due to emission of photons with frequency v. The number of 
electrons with energy E e , Q(E e ), generated per unit time in the GC due to 
XX-annihilation is given by 
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Q( Ee) =±pj:<* v > F ( d -^i) do) 

where x = E e /m x and dN e (x) / dx is denned in the same way as dN 1 (x)/dx. 

In the calculation of radio flux we assumed the magnetic field in the GC 
to be H = 10" 2 - 5 G (see below). 

We shall now briefly review the observational properties of the core of 
the GC following the references [12-13]. A violent activity is observed within 
a region of size ~ 1 pc, which can be powered by xx-annihilation. The to- 
tal bolometric luminosity, mostly due to the optical and UV radiation, is 
Lboi — 4 • 10 40 erg s _1 . The center of activity coincides with the com- 
pact non-thermal radio source Sgr A*. It has an unusual rising spectrum 
F v ~ l(i//10 GHz) 1 / 4 Jy in the frequency range 1 GHz < v < 86 GHz. Its 
luminosity is L v = 1.3 • 10 34 erg s -1 at wavelength A > 3 mm. The source is 
surrounded by a more extended radio halo. 

The X-ray radiation (0.5 — 4.5 KeV) from GC reveals [14] several com- 
pact sources one of which, with luminosity Lx — 1.5 • 10 35 erg s _1 , coincides 
with Sgr A*. The extended X-ray component around this source has lumi- 
nosity L x — 2.2 • 10 36 erg s -1 . 

High energy gamma-radiation from GC is not detected. The COS B 
upper limit at E 1 > 300 MeV is F 1 < 4-10 -7 cm -2 s _1 . There are indications 
of the presence of a powerful flux (F 1 ~ 10 -3 cm -2 s _1 ) in the form of 
0.511 MeV e + e~-annihilation line from GC. The core is submerged into a 
dust cloud, probably in the form of a ring with inner and outer radii of 
2 pc and 4 pc, respectively. The dust absorbs a considerable part of the 
radiation from the core and reradiates it in the form of IR radiation. The 
phenomenological model of the core is developed in refs. [12,13]. It has 
been shown that the observations provide an overdetermined set of relations 
for the properties of the core; these are satisfied by a unique and plausible 
solution [12]. In particular, the mean magnetic field in the core is found to 
be H ± = 10" 2±0 - 5 G. 

4. Results and Conclusions 

Let us turn now to the presentation and the discussion of our results. In 
Fig.l we report the values of the neutralino relic abundance in the form of a 
scatter plot obtained by varying Mi and \i in the ranges < Mi < 6 TeV 
and —3 TeV < \x < 3 TeV, respectively. The other parameters are taken 
at the following values: tan/3 = 8, rrih = 50 GeV, m? = m = 3 TeV (m? 
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denotes the mass of the sfermions, taken here to be degenerate). We notice 
that in this representative point of the parameter space many neutralino 
configurations provide large values for Q x h 2 ; for many configurations the 
relevant Vt x h 2 value is even above the cosmo logical bound Vth 2 = 1. This 
situation occurs for gaugino-dominated compositions, since in this case x~X 
annihilation proceeds mainly through /-exchange, but this process is now 
strongly hindered due to the large value assigned to mj. 

In Fig. 2 a scatter plot is displayed for the integral gamma-ray flux F 7 (> 
300 MeV) evaluated from eq.(8). We remind that in the calculation of F 7 
the rescaling for p l x c = ^xPdm ^ s use d- Fig. 2 shows that a large fraction 
of ^-configurations with m x < 100 GeV are excluded by the COS B upper 
limit F 7 (> 300 MeV) < 4 • 10" 7 cm" 2 s" 1 [11]. For the configurations 
allowed by this limit, the graph does not show the values of VL x h 2 , i.e. it 
does not specify whether the neutralino can provide Ocdm ~ 0.7. This 
informations is added in Fig. 3, where in a /j - M2 diagram are displayed 
the x-compositions which satisfy the two following requirements: a) Vt x is 
within the range: 0.5 < fi x < 1 which, taking h = 0.7, implies 0.25 < 
£l x h 2 < 0.5; b) the relevant radio flux, as evaluated from eq.(9), does not 
exceed the experimental upper bound F u < 1 • 10 -23 erg cm -2 s _1 Hz at 
v = 1 GHz [12,13]. The neutralino configurations that meet these properties 
are denoted in Fig. 3 by squares: they correspond to gaugino-dominated 
compositions with masses in the range 100 GeV < m x < 500 GeV. For 
these ^-configurations the relevant gamma-ray flux is 5 • 10~ 9 cm -2 s -1 < 
F 1 < 2 • 10~ 7 cm -2 s _1 , i.e. this source is detectable with existing or future 
devices. 

Fig. 4 displays the scatter plot of F 1 vs O x restricted to those configu- 
rations with mass 20 GeV < m x < 100 GeV, which satisfy the observational 
upper bound on radio-emission. The values fi x < 0.1 imply that CDM 
is composed mostly by some other particle(s), with neutralinos providing a 
small part of the CDM density. Fig. 4 shows that in this case the gamma-ray 
flux is marginally detectable down to fi x ~ 3 • 10~ 4 . It means that the NGS 
model can be tested with the help of ^-annihilation even if the neutralino 
constitutes only a small part of CDM. 

The general features of the previous discussion remain valid when we 
assign to the free mass parameters, mj rrih, values different from the pre- 
vious ones. For instance we could assign to mj the smallest possible value 
consistent with the LEP lower bound > 45 GeV [7] and with the con- 
dition that the neutralino in the lightest SuSy particle, and take for in- 
stance = 1.2 m x , when m x > 45 GeV and = 45 GeV otherwise. 
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In this case lighter sfermions make the x~ X annihilation for gaugino con- 
figurations much more efficient than before. As a consequence, now the 
fluxes for gamma-radiation and for radio-emission, which are proportional 
to < (tv >o /(< civ >int) 2 are smaller than in the previous case (where 
sfermions were assumed to be very massive). However, our calculations indi- 
cate that a situation similar to that depicted in Fig. 4 still persists, i.e. in a F 1 
vs fi/i 2 plot many neutralino configurations densely populate a region where 
O ~ 0.01 — 0.1 with values of F 7 only very slightly below the experimental 
upper bound. 

As a conclusion, we claim that the density distribution (1) of DM in 
our Galaxy results in a pointlike source of high-energy (E 1 > 100 MeV) 
gamma-radiation with fluxes of order ~ 10~ 7 — 10~ 8 cm -2 s -1 , detectable 
with the help of present and futures techniques. 
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Figure Captions 



Figure 1: Scatter plot of the neutralino relic abundance as a function of 
the neutralino mass. This scatter plot has been obtained by varying M 2 and 
H in the ranges < M 2 < 6 TeV and -3 TeV < fi < 3 TeV respectively. The 
other parameters are taken at the following values: tan/3 = 8, rrih = 50 GeV, 
rrij = m = 3 TeV {mj denotes the mass of the sfermions, taken here to be 
degenerate) 

Figure 2: Scatter plot of the integral gamma-ray flux F 7 (> 300 MeV) 
as a function of the neutralino mass. The ranges for M 2 and ji are the same 
as in Fig.l. Also the other parameters are taken at the same values as in 
Fig.l. 

Figure 3: ji— M 2 diagram where the squares denote the x-compositions 
which satisfy the two following requirements: a) fi x is within the range: 0.5 < 
fi x < 1 which, taking h = 0.7, implies 0.25 < £l x h 2 < 0.5; b) the relevant 
radio flux, as evaluated from eq.(9), does not exceed the experimental upper 
bound F v < 1-10 -23 erg cm -2 s _1 Hz at v = 1 GHz [12,13]. The gamma-ray 
flux for these configurations is 5 • 10~ 9 cm -2 s _1 < F~, < 2 • 10~ 7 cm~ 2 s _1 

Figure 4: Fig. Scatter plot of F 1 vs fi x restricted to those configura- 
tions with mass 20 GeV < m x < 100 GeV, which satisfy the observational 
upper bound on radio-emission. 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402215vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402215vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402215vl 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402215vl 



